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p Abstract 

General relativistic effects have been shown to increase the energy deposition rate 
due to the process vv — > e + e~ in supernovae and neutron stars. In this paper we 
study the effect of inclusion of the rotation of the star in the general relativistic 
treatment. We show that inclusion of rotation results in a reduction in the heating 
\ rate as compared to the no rotation case. 
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Oh 1 Introduction 

6 

The neutrino pair annilation reaction (vv — > e + e ) is one of the important pro- 
cesses in understanding the energy transfer from a hot proto neutron star to the 
outer layers of a supernova [|TJ. It was pointed out in 0, |3| that this process can 
give a great boost to the delayed neutrino heating mechanism of Bethe and Wil- 
son j|] in addition to the i/-nucleon capture reactions. The annihilation process 
is useful as this can continue giving energy to the "radiation bubble" till there is 
an emitted neutrino flux ||. Apart from type II supernovae this process is also 
important for collapsing neutron stars ||, close neutron star binaries in their last 
stable orbit ||. In particular, this process has been considered to be one of the 
possible sources of energy for gamma-ray bursts @. The reaction efficiency of this 
process was calculated in ^ ]j| using Newtonian gravity. However strong grav- 
itational effect in supernova and collapsing neutron star envioronments render a 



general relativistic calculation necessary [jTO|, [II], [12], [T^, |TJJ]. It has been shown in 
|12| that the energy deposition rate increases by about a factor of 4 in supernovae 
and by a factor of 30 in neutron stars if general relativistic effects are taken into 
account. In this paper we include the effect of rotation in the general relativistic 
treatment of [|I^] and compare how much this changes the reaction efficiency of 



the neutrino pair annhilation process over the case discussed in |12| . 



1 



2 Geodesies in the External Field of a Slowly 
Rotating Object 

The metric outside a slowly rotating star as given by the approximate Hartle- 
Thorne solution, with only the dipole corrections on a static star is given by |15| 

ds 2 = _(l_2s+ 2 4)dt 2 + (l - 2 -^ + ^Y 1 dr 2 + r 2 d9 2 

V r r (1) 

+r 2 sin 2 9 (d<j> - %dt) 

where M is the mass and J the specific angular momentum. Here and through out 
the paper (unless otherwise mentioned) we use the geometrised units (G=c=l). 
Using the Lagrangian appropriate to the motion in the equatorial plane (6 = 7r/2) 



2C -(l-^-^)t 2 +(l-^ + ^)^r 



(2) 



the generalised momenta may be written as 

Pt = -(i-*?-¥)i-¥4> = -E 



P<t> = -^ + r 2 = L (3) 



Pr = (l-^ + ^) _1 r 
and thereby the Hamiltonian 

2M 2J 2 



-l 



2H = -Ei + L<j) + I 1 - + r 2 = 5x (4) 

wherein the constant 5x — 1 for time-like geodesies and 8\ — for null geodesies 

Solving for t and (ft, one gets 

U* = i = (E- 2 £)(l- 2 f + 2 £ 



16]. 



JJ<t> — 1 - \(l _ 2M _ 2J 2 \ L _|_ 2J£l (l _ 2M_ , 2J 2 \ 1 



(5) 



As we will be interested in zero rest mass particles we consider the null geodesies 
5i = 0, and thereby obtain r from (4) to get the equation 



which in the limit J — > reduces to the well-known Schwarzschild form [17 

I I 2 A. I ^ 



1 dr 



iy 2 
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where L/E — > b, the impact parameter for a massless particle. 



Introducing the Local Lorentz tetrad X^ a \ as given by 
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(upper index refers to rows and the lower to columns), one can define the angle 
r between the trajectory and the tangent vector in terms of local radial and 
longitudinal velocities 



tan6 T 



v 3 



I ^_ 2M_ , 1j' N 



(9) 
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As the local velocity = U^/U 1 , using (5), one can get 



dr N 



2M 2J 
1 + — 



2\ 3/2 r 



2M 2J 2 



1 2J 

r 6r 2 



n -1 



tan Or 



(10) 



Eliminating f^j between (6) and (10) and simplifying, one gets the impact pa- 
rameter b to be 

1 -2M/r + 2J 2 /r 4 ) 1/2 ' ~' 



2J 



r cos r 



11) 



For a neutrino emitted tangentially (On = 0) from a neutrinosphere at R 



which reduces to 



cos 9 r 



R 



1 - 2M/r 
\ 1 - 2M/R 



for J = ||12|| . The minimum photosphere radius for the non-rotating case was at 
R = 3M. For the present case the minimum radius and the corresponding impact 
parameter may be found from the roots of the system of equations 



R A - (3 - 6 J/6) R 3 - 6 J 2 



and 







R b - 6"i? 4 + 26(6 - 2J)R 6 - 26 2 J 2 







(13) 



(14) 



obtained from the effective potential for the particle in circular orbit. In the above 
equations we have expressed all the quantitites in units of stellar mass M. The 
maximum value of the rotation parameter that we take is J/M 2 = 0.3 as the 
metric used is valid for slowly rotating configurations. It may indeed be verified 
that this value J/M 2 (= 0.3) is about the critical value for the case of a millisecond 
pulsar PS| . In Table 1 we give the values of neutrinosphere radius (R u ) (~ last 
photon orbit) and the impact parameter (6) for different values of J/M 2 by solving 
simultaneously the eqs. ( |T3"D and ([14]). For every value of J/M 2 there exists atleast 
one real root which is always located at a distance < 3M, giving the neutrinosphere 
radius to be less than that in the case of Schwarzschild solution (the non-rotating 
case). In Table 1 we also give the corresponding event horizons (Reh), obtained by 
solving the equation g u = 0, for various values of J/M 2 . In the following section 
we calculate the energy deposition rates at a distance r > R v . 



Table 1: The values of the event horizon (Reh), neutrinosphere radius (R u ) and 
the impact parameter (b) for different values of the rotation parameter J/M 2 . 



J/M 2 


Reh/M 


Rv/M 


b/M 


0.1 


2.00249 


2.8771 


4.9855 


0.2 


2.00985 


2.7355 


4.7483 


0.3 


2.02178 


2.56546 


4.4715 



3 Energy Deposition Rates 

The energy depostion rate/unit time/unit volume by this process is given in general 
by§ 

i= 1D ffj^\ mr))Mr) (15) 
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where Gp is the Fermi coupling constant, and D = 1 ± 4sin 2 9 w + %sin A 9w with 
sin 2 8 w = 0.23 and the + sign is for electron type neutrinos and antineutrinos and 
the — sign is for the muon and tau type neutrinos and antineutrinos. T(r) is the 
temperature measured by the local observer and 0(r) is the angular integration 
factor. A general relativistic treatment requires the incorporation of gravitational 
red shift in T(r) and the effect of bending of the path of neutrinos in 0(r). In 
terms of the unit direction vector fl v and the solid angle subtended dQ v , 0(r) can 
be written as 



9(r) 



1 



Q, v .Q, 9 ) 2 dQ, u dQ,p 
47T 2 / / [1 - 2/i^/ip + nlnl + -i 

Jx Jx Z 



(16) 

//„ 2 )(1 - i4)]d^ud^ m (17) 



where \i = sin9, Q = (//, a/1 — fi' 2 cos(j), a/1 — fi' 2 sin(f)) and dVL = cos9d9d(p. This 
simplifies to 

2tt 2 

0( r ) = _(l-x) 4 (a; 2 + 4a; + 5) (18) 



where x = sin9 r . With 9 r defined in eq. 

#y4(l- ^. + 1 



of the previous section one gets 

1/2 
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(19) 



This reduces to the expression given in for J = 0. 

The neutrino temperature varies linearly with redshift and T(r) in eq.(|l5D at a 
radius r is related to the neutrino temperature at the neutrinosphere radius R as 



T(r) 



i 2M 2J 2 

\ 1 _ 2M _ 2J£ 
\ 1 r r 4 



T(R) 



(20) 



The total amont of local energy deposited by vv — > e + e~ for a single neutrino 
flavour for a rotating star can be defined as (in analogy to the non-rotating case 
discussed in [ |12]] ) 

47rr 2 dr 



Q 



R 



2M 



2J2 
^4 



(21) 



where g is defined in eq.(|15|) with T(r) and 0(r) defined in eqs. (PP]) and (Tl6| ) 
respectively. This can be simplified to obtain 



Qr.j 1.09 x K)-'^ (^,^ JDI,! 



9/4 p-3/2 



f22l 
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where Q51 and L 51 (luminosity) are in units of 10 51 ergs/sec, Rq is the radius in 
units of 10 km and 

9/4 



V dy 



( M 


R 2 ) 


= 3(l-¥ 


2J 2 \ 


\ R ' 






iJ 4 y 



r (1- l) 4 (x 2 + 4x + 5) • -7- 

1 C 1 - 



2M 2 J 2 
r r 4 



where ?/ = r/i2. For J — > this reduces to .F(^r) of |12| and the Newtonian limit 
is obtained by taking M — >■ H. 



4 Results and Discussions 

In fig. la and fig. lb we plot the ratio Q(%, ^)/<?Newt = if> if) vs R/M. 
In fig. la we plot the ratio from R/M = 5 to R/M = 10 which is relevant for 
a type II supernova while in fig. lb we plot the ratio at smaller values of R/M 
which is relevant for collapsing neutron stars. The J —>■ curves correspond to the 
case without rotation for which the general relativistic effects enhances the energy 



deposition rate by a factor of 4 at R = 5M |L2| and almost by a factor of 30 for 
R = 3M. As one includes the rotation of the star there is a drop in the energy 
deposition rate. The reduction is more for smaller values of R/M and higher values 
of J/M 2 . If we take J/M 2 = 0.3 then at R/M = 3 we get a reduction by « 38 
% whereas for R/M = 5 the reduction due to rotation is ~ 9%. With rotation 



the apparent angular size of the star seen by the neutrino decreases (see eq. dl2|) 
thus decreasing the probability of head on collision. This results in a drop in the 
heating rate. 

In fig. 2 we plot dQ/dr vs the radius for three values of R/M (5,7 and 10). For 
each case we plot dQ/dr for various values of the rotation parameter J/M 2 . We 
also plot the Newtonian case for which M=0 and J=0. The heating rate is seen 
to fall sharply as the radius increases. The inclusion of general relativistic effects 
enhances the heating rate but introduction of rotation reduces this. The reduction 
is more pronounced for smaller values of R/M and higher values of J/M 2 . The 
resultant heating rate after including the rotation is however still higher than the 
Newtonian value. 

To conclude, in this paper we have extended the general relativistic calcula- 
tions of the neutrino heating rates due to neutrino pair annihilation to include the 
rotation of the star. We find that the effect of including rotation is to reduce the 
heating rate over the no rotation case and the reduction can be as large as 38%. 
The effect is more pronounced for smaller values of the ratio R/M, in the range 
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relevant for collapsing neutron stars. 
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Figure Caption 

Fig. I. The ratio of the general relativistic energy deposition rate to the New- 
tonian rate for various values of the rotation parameter. The solid line is for no 
rotation, the big dashed line is for J/M 2 = 0.1, the small dashed line is for J/M 2 
= 0.2 and the dotted line is for J/M 2 = 0.3. 

Fig. 2. dQ/dr as a function of radius for different values of M/R and J/M 2 . 
Also shown is the variation for the Newtonian case (M=0,J=0). 
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